An Improved Description of the Dielectric Breakdown in Oxides Based on a 

Generalized Weibull distribution 
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In this work, we address modal parameter fluctuations in statistical distributions describing 
charge-to-breakdown (Qbd) and/or time-to- breakdown {tBD*) during the dielectric breakdown 
regime of ultra-thin oxides, which are of high interest for the advancement of electronic technol- 
ogy. We reobtain a generalized Weibull distribution (q- Weibull), which properly describes (tsc) 
data when oxide thickness fluctuations are present, in order to improve reliability assessment of 
ultra-thin oxides by time-to-breakdown (tso) extrapolation and area scaling. The incorporation 
of fluctuations allows a physical interpretation of the q- Weibull distribution in connection with 
the Tsallis statistics. In support to our results, we analyze tso data of Si02-based MOS devices 
obtained experimentally and theoretically through a percolation model, demonstrating an advanta- 
geous description of the dielectric breakdown by the q- Weibull distribution. 

PACS numbers: 71.20.Nr,71.30.+h, 77.22.Jp, 85.40.-e,02.50.Cw 



I. INTRODUCTION 

The International Technology Roadmap for Semicon- 
ductors indicates the need to decrease the SiC>2 gate 
thickness to less than 3 nm to pursue the continuity of 
the Moore's Law Q. Research efforts to maintain this 
trend have been focusing on: (?) the advantageous sub- 
stitution of SiC>2 by high dielectric constant materials to 
allow equivalent gate dielectric oxide thickness bellow 1 
nm 0; (ii) a better understanding of the fractal regime 
of ultra-thin gate dielectric oxide growth, where spatial 
inhomogeneities or fluctuations are important 0, Q con- 
siderations for the control of dielectric breakdown condi- 
tions. Consequently, the reliability of ultra-thin oxides 
for ultra-large scale integration is one of the most impor- 
tant concerns in the domain of electronics miniaturiza- 
tion nowadays 0, IE • 

Two quantities which are experimentally measured 
to assess the reliability of a metal-oxide-semiconductor 
(MOS) device are: the charge-to-breakdown Qbd, de- 
fined as the time-integrated current density which flows 
until breakdown occurs, and the time-to-breakdown 
tsD, defined as the samples' lifetime. These quan- 
tities are statistically distributed and are usually as- 
sumed to be Weibull distributed In particular, 
accurate assessment of reliability of ultra-thin oxides 
is seriously affected by fluctuations of oxide thickness. 
These fluctuations cause a deviation from the Weibull 
distribution, which is easily observed by a curvature 
in a graph of ln[— ln[l — -F(i)]] versus ln[t], being 



F{t) = f Q P(x)dx the cumulative distribution. Extrap- 
olation of the slope of this plot from higher failure per- 
centiles to lower failure percentiles can lead to serious 
errors in reliability assessment, as shown by Wu et al. 



Deviations from a given distribution can be generated 
by taking scale mixtures of it |9(. For instance, the scale 
parameter of the exponential and Gaussian distributions 
can be averaged by a gamma distribution. In these two 
cases, the deviation induced by the fluctuation in this pa- 
rameter has been successfully employed in the description 
of physical situations, such as those related with nonex- 
ponential decay and turbulence [T3.[T3] |. Further- 

more, this average process has also been connected with 
the entropic parameter q in the context of Tsallis statis- 
tics [jjj, [HI lQ • By taking the fluctuations of the modal 
parameter of the Weibull distribution into account, we 
obtain a generalized Weibull distribution which properly 
describes tsD data when oxide fluctuations are present, 
thus improving reliability assessment of ultra-thin oxides 
by tsD extrapolation and area scaling. This generalized 
Weibull distribution, without focusing fluctuations, was 
discussed and applied in Ref. 0] to other systems. 



In the following section we present the generalized 
Weibull distribution employed here. The application of 
this distribution to the description of the dielectric break- 
down in oxides is considered in section 3. The last section 
is dedicated to our conclusions. 
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FIG. 1: (a) Graph of ln[-ln[l - F(t)]] versus t BD for 4000 
ultra-thin oxide MOS capacitors. The data were obtained 
from the work of Wu et al. 5]; (b) the same data as in 
(a) but depicted in a g- Weibull plot in which a correct tBD 
extrapolation is very clear. 



II. g- WEIBULL DISTRIBUTION 

The Weibull distribution often used in the statistical 
description in the study of the time-to-breakdown in elec- 
tronic devices is given by 



Pi(x) = 



x \x 



0-1 



exp 



(1) 



where the parameter xq is the modal value that corre- 
sponds to 63, 2% of the samples' lifetime. The parameter 
(3 is the Weibull slope and defines the time-to-breakdown 
spreading degree. In order to incorporate the referred 
fluctuation, we consider a Weibull compound distribu- 
tion p(x; 9) — f39x l3 ~ 1 exp(— 9x@). Here, 9 — x^ 13 , with 
(3 fixed and a:o changing so that the new parameter 9 
is a stochastic variable ruled by the gamma probability 
density function 



p(9) 



Sr r-l e -80 



r(r) 



?,r,<5>0). (2) 



FIG. 2: (a) Graph of ln[-ln„[l - F d (t)]] versus ln(t) for a 
6.9 nm sample oxides subjected to a Eox = 12.7 MV/ cm ap- 
plied field. The samples have areas of 10 -4 and 10~ 5 mm 2 . 
The data were obtained from the work of Teramoto et al.) 
pill ]: (b) the same data as in (a) but depicted in a q- Weibull 
plot showing improved area scaling. 



Therefore, the average in the stochastic parameter 9, 

/•oc 

{p{x;9)) e = / p{x;6)p{6) d0 , (3) 
Jo 

leads to 

(p(x; 9)) e = r5 r px - x {x + S)-^ r+1 \ (4) 

Note that by naming this average P q i(x), it can be 
written as 



x \xoJ 



(5) 



where we used the g-exponential function denned as 
ex P(? , (-a a;) = [l-(l-g') ai] 1 /! 1 -''' if l-(l-g') ax > 0, 
and exp 9 /(— ax) = if 1 — (1 — q') ax < 0, with 

q' = (2 + r)/(l + r) and x% = [5/(r + 1)]. 

By comparing Eq. (JIJ with Eq. JSJ) we verify that P q i (x) 
is a generalization of the Weibull distribution in the same 
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Sq = (1 — 2j=i Pi )/(! i s maximized subjected to 
appropriate constraints [jj]]- Here, pt is the probability 
of the ith state, W is the number of accessible states, 
and q is a real parameter that rules the degree of gener- 
alization of the theory (when q — > 1 we recover the usual 
entropy). On the other hand, the g-exponential func- 
tion can be obtained from the average of the exponential 
function by using the gamma distribution, Eq. J2J) 18] . 
Thus, the parameter q occurring in the Tsallis statistics 
is shown to be entirely induced by the fluctuations of the 
parameter of the usual exponential distribution [lOj . Fol- 
lowing this interpretation, the average process employed 
to obtain P q i (x) characterizes an implicit connection with 
the Tsallis statistics. 
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FIG. 3: (a) Graph of ln[-ln[l - F d (t)]] versus ln(t) from 
the percolation simulations as performed by Sombra et al. 
PH for gate (squares) and substrate (dots) injection adjusted 
with Weibull distributions; (b) the same data as in (a) but 
depicted in a q- Weibull plot showing improved fitting. 



sense that the g-exponential function generalizes the ex- 
ponential one. In view of this observation and following 
Ref. Q , we refer to P q < (x) as q- Weibull distribution. We 
also remark that P q > (x) contains a broad class of distri- 
butions as particular cases. In fact, in the limit q' — > 1 it 
reduces to the Weibull distribution, for j3 — ► 1 it gives the 
q-exponential distribution, and when q 1 — > 1 and j3 — > 1 
it becomes the exponential distribution. Furthermore, 
from Eq.(f4J| we verify two power law regimes, one for 



small 



0-1 



and another for large x, x 



Notice 



also that L P q i(x)dx exists and is equal to one when 
q' < 2. 

Before addressing our discussion to the application of 
the q- Weibull distribution concerning experimental and 
simulation results, we give a connection of this distri- 
bution with the Tsallis statistics. The first step to- 
wards verifying this connection is based on the fact that 
the g-exponential can be viewed as a signature of the 
Tsallis statistics, since it basically replaces the exponen- 
tial in the canonical distribution. In fact, this canon- 
ical distribution is obtained when the Tsallis entropy 



III. APPLICATION TO DIELECTRIC 
BREAKDOWN 

To apply our distribution, Eq. (JSJ), to investigate ex- 
perimental and numerical simulated data, we consider 
the cumulative distribution for P q >(x), i.e., 



F q (t) = 1 - ex P(J 



(6) 



with q = 1/(2 — q') and a = xq/(2 — q') 1 ^ . In connection 
with the ^-exponential function, the (/-logarithm function 
is usually defined as 



A-q _ i 



thus \n q [exp q (x)] = exjp q [ln q (x)] ■ 
we obtain that the graph of ln[- 
ln(£) gives a straight line since 



(7) 
q 

- x. By using Eq. 
ln g (l — F q (t))] versus 



ln\ 



-ln q [l - F q (t)]] =(3ln(t) - pin(a). 



(8) 



Therefore, if a cumulative distribution relative to the di- 
electric breakdown for a given set of data, F d (t), is well 
described by the q- Weibull distribution, then the q pa- 
rameter present in the q-logarithm is obtained through a 
linear adjustment, and the parameter /3 is the slope of the 
graph. This procedure introduced here will be referred 
as a q- Weibull plot. 

To compare the usefulness of the Weibull and the q- 
Weibull distributions, we analyze experimental |19| and 
simulation results. Figure ^a), which is adapted from 
the work of Wu et al. j5| , shows a curvature in a graph 
ln[— ln[l — versus \n[t]. The variation of this slope 

is associated with oxide thickness fluctuations. In fact, 
Wu et al. were able to generate this slope variation 
by performing simulations in which the oxide thickness 
obeyed a Gaussian distribution. Extrapolation of the 
slope of this plot from higher failure percentiles (conven- 
tional experimental window) to lower failure percentiles 
leads to an error in tsD projection and, consequently, 
in projection of reliability. It is necessary to use large 
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and time-consuming sample sizes (thousands of devices) 
in order to obtain tsD distributions down to lower per- 
centiles and avoid erroneous projections. The same tsD 
data can be correctly fitted in the entire range of per- 
centiles by the g- Weibull distribution, as shown in Fig. 
Gib). In the q- Weibull plot, the data follow a straight 
line. We conclude that, by taking the fluctuations of 
the modal parameter of the Weibull distribution into ac- 
count, one obtains a statistical distribution which com- 
pletely describes tsD data when oxide fluctuations are 
present. Therefore, the q- Weibull distribution allows im- 
proved tBD projection for ultra-thin oxides and does not 
require a very large number of samples, i.e., correct tsD 
extrapolation can be performed in the conventional ex- 
perimental window. This is a striking benefit brought by 
the q- Weibull distribution. 

In addition to correct tsD projection, we have found 
that the q- Weibull also improves area scaling, another 
important consideration for reliability studies |15|. In 
Fig. |2I a), we have used the data from the work of Ter- 
amoto et al. [20|. They studied the time-dependent 
dielectric breakdown for 6.9 nm thin SiC>2 oxides with 
samples areas of 10 -4 and 10~ 5 mm 2 , and subjected to 
a Eox = 12.7 MV/ cm applied field. It is clear, when 
comparing Fig. Efa) and Fig. Efb), that the q- Weibull 
distribution gives a better area scaling than the Weibull 
one, another remarkable benefit. 

When tBD data are plotted in a limited range of per- 
centiles, the differences between the q- Weibull and the 
Weibull distributions become less apparent, but can still 
be seen. As an example, we consider the simulations 
performed by Sombra et al. 21]. They developed a per- 
colation model to describe the dielectric breakdown of a 
MOS capacitor investigating effects of bias polarity, ox- 



ide film thickness and electric field strength. The hot 
electron injection is either through the gate or the sub- 
strate. In Fig. EJa) we have performed a linear fit of the 
simulation data of Sombra [21| with the Weibull distri- 
bution, while in Fig. EJb) we depict the same data but 
with a linear fit through the q- Weibull distribution. One 
can see once more that the q- Weibull distribution gives a 
better fit than the Weibull one, for small and large tsD- 



IV. CONCLUSIONS 

In summary, we conclude that: (i) fluctuation on the 
modal value is a mechanism to explain deviations from 
the Weibull distribution in reliability studies of electron- 
ics devices; (ii) this fluctuation leads to a generalization 
of the Weibull distribution (q- Weibull distribution) and a 
connection with the Tsallis statistics; and (Hi) the gener- 
alized Weibull distribution leads to a better adjustment 
of experimental and simulated data in comparison with 
the Weibull one, giving improved tsD extrapolation and 
area scaling. Finally, the results presented here can be 
very useful for an improved description of the dielectric 
breakdown in the high dielectric constant materials Q 
which will be present in future generations of advanced 
MOS devices. 
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